Measuring charge based quantum bits by a superconducting single-electron transistor 
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Single electron transistors have been proposed to be used as a read-out device for Cooper pair 
charge qubits. Here we show that a coupled superconducting transistor at a threshold voltage is 
much more effective in measuring the state of a qubit than a normal-metal transistor at the same 
voltage range. The effect of the superconducting gap is to completely block the current through 
the transistor when the qubit is in the logical state 1, compared to the mere diminishment of the 
current in the normal-metal case. The time evolution of the system is solved when the measuring 
device is driven out of equilibrium and the setting is analysed numerically for parameters accessible 
by lithographic aluminium structures. 
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Nanoscale devices such as Cooper pair boxes or cou- 
pled quantum dots have been suggested as scalable and 
integrable realisations of quantum bits. The two logical 
states of a qubit are the different charge states, or in the 
case of a flux qubit, flux states of the system. There 
are several proposals for quantum gates f 1 and for in- 
terqubit couplings as well as for measuring devices 
@j Hi Hi 0- Permanently coupled normal-metal transis- 
tors have been suggested as a device for measuring the 
state of a Cooper pair charge qubit @, |). Also, a su- 
perconducting SET in the Josephson current regime (low 
biasing voltage) has been experimentally tested flio|| . In 
this work we show that in the regime of high biasing volt- 
ages, the superconducting SET Jll], |l^| leads to a highly 
efficient quantum nondemolition measurement Jl^| due 
to the blocking effect of the gap. 

The setting is shown in Fig. [j]. In the upper part, the 
Cooper pair box forms the qubit, its state characterised 
by the number of excess Cooper pairs in the box, n. In 
the lower part, the superconducting single-electron tran- 
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Figure 1: A Cooper pair charge qubit is capacitively coupled 
to a measuring single-electron transistor. The quantum num- 
bers n, N and m are explained in the text. The SET is sym- 
metrically biased with voltage Vbias, and the voltage V sa t is 
the gate voltage. The energy scales are determined by the in- 
teraction capacitance Cmt, the tunnel junctions' capacitances 
Ct and the gate capacitance C se t- 



sistor is capacitively coupled to the qubit with its state 
characterised by the excess charge on the island, eA^. In 
addition, the quantum number to counts the number of 
charges passing through the SET in left-to-right direc- 
tion. Without a biasing voltage Vbias across the SET, 
there is no dissipative current and no information is re- 
ceived. Moreover, in the absence of a dissipative environ- 
ment, no dephasing of the composite system will occur 
and quantum operations on the qubit can be performed. 

In order to perform the measurement, a biasing volt- 
age is applied. As different qubit eigenstates corre- 
spond to different conductance in the SET, by observing 
the current one receives information on the state of the 
qubit. The time needed for the current to give the es- 
sential information is called the measurement time. The 
back action caused by the SET dephases the qubit and 
eventually destroys also the logical state of the qubit. 
The corresponding time scales are called the dephas- 
ing and the mixing time, respectively. For a good non- 
demolition measurement of the logical state (|a| 2 and \b\ 2 
in a|0) + b\l)), one expects to have a much longer mixing 
than a measurement time scale. We show that for a su- 
perconducting SET, in a parameter range accessible by 
aluminium structures, the ratio between the mixing and 
measurement time scales is in excess of 600. 

The total Hamiltonian consists of three parts: the 
Hamiltonians of the SET, the qubit and the interaction 
H sct , -f/qb and -ffmt, respectively. The SET Hamiltonian 
is defined as 



E set (N -Qsetf + H T + H L +H R + H 1 , (1) 



where E sct is the charging energy, Q se t the gate charge 
of the transistor and 
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(2) 

describes quasiparticle tunnelling within the SET. The 
phase differences cf> and \& are the conjugate variables of 
N and to, respectively. Thus, the operator e 1 ^ (e 1 *) cor- 
responds to quasiparticle tunnelling onto the island (right 
electrode) of the SET, increasing the quantum number N 
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m) by one. The last three terms in Eq. (|l]) are defined 



as H r = J2k<r e fc"-4a c L> (r = L, R, I), where a labels the 
transverse channels including the spin and k labels the 
wave vector within one channel. These describe the non- 
interacting electrons in the left electrode, the right elec- 
trode and the island, respectively. Cooper pair tunnelling 
is excluded from the SET Hamiltonian due to the high 
biasing voltage. However, this approximation is not valid 
for a low voltage, and the contribution from the Cooper 
pair tunnelling should be added when analysing the effect 
of the SET on the qubit during logical operations. 
The qubit Hamiltonian is defined as 
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(3) 



where E q ^ is the charging energy, Q q b qubit 's gate charge 
and the last term describes the transfer of Cooper pairs 
to and from the box (e ie |n) = \n + 1}). The quasiparti- 
cle tunnelling in the qubit is suppressed by the Coulomb 
blockade, and the microscopic degrees of freedom have al- 
ready been integrated out. The interaction Hamiltonian 
describes the Coulomb interaction between the qubit and 
the transistor and is defined as 



H, 



E- m ,nN, 



(4) 



where E- lnt is the charging energy. All charging energies 
E mt , E sot and -E q b are determined by the capacitances of 
the system, and the gate charges Q q b = —eV n /E q b and 
Qsot = — eVjv/2-E l set depend on the effective gate voltages 
V n and Vn that, for a symmetric bias, are determined by 
the gate voltages V q h and V sct as given in M . 

Concentrating on the values of Q q b around the degen- 
eracy point Q q b = |, only the low-energy charge states 
ri = and n = 1 are relevant. This assumption is re- 
quired for the Cooper pair box to constitute a quantum 
bit. In the qubit's two-state approximation, the diago- 
nalised operator becomes H q h — —\AEa z , where 



AE= J[E qh (l-2Q 



qb) 



(5) 



while the operator n becomes nondiagonal n = 
h (1 — cos(ry)cr 2 — sm(r))a x ), where the mixing angle is 
given by tanry = Ej/ [-Eqb(l — 2Q q b)]. Assuming the 
gate charge Q q b sufficiently different from the degeneracy 
point 0.5, the mixing angle ry is small in the charge-qubit 
approximation Ej <C -E q b- By rearranging the operators, 
the final form for the total Hamiltonian can be written 
as H = Hq + Hi , where 

H = H L + H R + H 1 + H qh + E sct (N - Q sct f + H int . (6) 

We analyse the measurement process by master equa- 
tion techniques. The master equation for the system 
reads 



doit) i . , ., „ 
^ + -[tfo,a(t)]=Tr w 



E(* - t')a(t')dt' ', (7) 



where the trace is taken over the microscopic degrees 
of freedom of the transistor's left and right electrodes 
and the island. The elements of the transition matrix 
E(t - t') = {[V, U{t - t') [V, •] U(f - t)]) are calcu- 
lated by using the diagrammatic technique developed in 
Q and ||. However, the rules for converting diagrams 
into integral equations must be changed to include the 
superconducting density of states 



\x\ 



N(x) = V* 



No 



if |x| > A 
otherwise, 



(8) 



where A is the superconducting gap. In zero-temperature 
limit, the matrix elements H(t — t') can be written as 



Ej, j (t - t') = F(t - i ')2 e -is(t-t')/fi ! 



(9) 



where F(t) — / °° N(x)e~^ +lt ' >x dx, and e is a cut-off 
exponent which ensures the convergence of the integral. 
With the superconducting density of states (||) , the func- 
tion F(t) is a so-called Basset function of first order |l(| 
which must be analysed numerically. 

By performing the Laplace transform on the master 
equation (0), the right-hand side becomes H(s)a(s). As- 
suming the density matrix a to change slowly in a time 
scale of H/E, the calculations can be restricted to the 
regime s <C E. In the normal-metal case, E(s) varies only 
slowly as a function of small s, and therefore the zeroth- 
order approximation is reasonably good. However, when 
the energy E/2 is close to the gap energy A, E(s) has 
a strong s-dependence which is approximately linear for 
small s, as shown by Fig. ||. Using the linear approxima- 
tion E(s) = a + b ■ s and performing the inverse Laplace 
transform, the right-hand side becomes (V + Ad/dt) a(t). 
Moving all the derivatives of the density matrix to the 
left-hand side gives as a master equation 



(1-A) 



da(t) 
dt 



h 



[a(t),H ]=Ta(t), 



(10) 



where F and A are tridiagonal matrices consisting of the 
zcroth- and first-order terms of the E(s) matrix, respec- 
tively. To guarantee the existence of the inverse of (1 — A) , 
the elements of the A coefficient are required to be small. 
This requirement is fulfilled when the tunnelling rate 
within the SET is small (that is, |T L |, |T R | < 1). Multi- 
plying equation ( [l0| ) from the left by (1 — A) -1 the final 
form of the master equation is obtained as 



da(t) 
dt 



(1-A)- 



[a(t),H ]+Ta(t) 



Ta(t). 



(11) 

This describes a Markovian evolution of the density ma- 
trix and a sequential tunnelling approximation, which is 
a reasonable assumption when the current in the SET is 
low. 
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Figure 2: The complex transition coefficient E(s + iE) plotted 
as a function of s with fixed E — 10.0 for different values of 
the gap A. The real components are drawn in solid curves 
and the imaginary parts in dotted ones. The corresponding 
gap values are in the order from the topmost curve to the 
lowest curve: A = 0.0, A = 3.0, A = 4.5, A = 4.9, A = 5.1 
and A = 5.3. For small values of s, the real components are 
nearly linear. 



The elements of the S(s) matrix are analysed by ex- 
plicitly writing the Laplace transform 



S( S ) 



/ Z(t)e- St dt= / F(t) 2 e { ~ s+lE)t dt. (12) 
Jo Jo 



The linear approximation is done by calculating S(s) for 
two values of s — 0.0 and s = 0.1 and fitting a line. 

Elements of the matrix T describing the transitions 
within the qubit are proportional to a small mixing an- 
gle e defined by tane = E- lnt smr]/(AE + E- mt cos?y). By 
approximating these elements by zeroes, the matrix T 
separates into four parts: one part describing the system 
when the qubit is in (diagonal) state 00, one part when 
the qubit is in state 11 and two parts for the nondiago- 
nal qubit states 01 and 10. The first two parts give raise 
to the measurement time as they describe the transis- 
tor's evolution for different qubit eigenstates. The co- 
herence of the qubit is described by the strength of the 
non-diagonal elements, and thus the rate at which the 
non-diagonal elements vanish gives the dephasing time. 
For the mixing time, one has to include the small mixing 
elements proportional to e. To ease the treatment be- 
low, a Fourier transformation is performed in m-space: 

The measurement time is determined by the first two 
parts of the matrix T. The equation corresponding to the 
qubit state 00 is 



d 
dt 



_N,k 
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N+l,k 

T 00 
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and the other part is a similar matrix equation for the 
qubit state 1 1 . These equations describe slowly damping 



(due to shot noise) conductance peaks that propagate 
in time. The measurement has been performed once the 
peaks can be distinguished from each other, and following 
the treatment in Jl| the measurement time t ms is defined 
as the time when the width of the peaks is smaller than 
the distance between their centers yielding 



(14) 



(Fi 



where the group velocities are Ti := Tl u 
T Ru ), with i = 0,1. The factors /* 
^r-)/(^l u + r_RiJ 2 are known as Fano factors [|3j. 

The dephasing time is determined by the evolution of 
the nondiagonal qubit states given by the remaining two 
parts of the matrix. Choosing k — has the effect of 
tracing out the ^.-variables and simplifying the corre- 
sponding submatrices 



dt Vcr, 
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(15) 



Again, these equations have plane wave solutions, and the 
corresponding eigenvalues determine the dephasing time 
of the system. The time corresponding to each eigenvalue 
A is defined as l/Re(A), and the dephasing time is chosen 
to be the longer of these two times. The shorter time 
describes partial dephasing. In the superconducting case 
these must be evaluated numerically. 

In order to analyse the mixing time of the system, also 
the mixing coefficients must be included in the T ma- 
trix. Once again, the SET degrees of freedom are irrele- 
vant and k is chosen to be 0. Now the matrix has eight 
eigenvalues: four of the eigenvalues describe the dephas- 
ing, two describe the measurement and one eigenvalue 
is zero (describing the trace preserving symmetry). The 
remaining real eigenvalue A m i x gives the mixing time as 
the t m i X := 1/A m ix, which must be analysed numerically. 

If the mixing time is very large compared to the mea- 
surement time, the measuring device disturbs the proba- 
bility amplitudes of the qubit (|a| 2 and \b\ 2 in a|0) + 
only a little. The uncertainties of the charge AQ and 
its conjugate variable (phase or flux) A<& are linked to 
each other by the uncertainty principle AQ • A$ > H/2. 
According to this principle, if the precision of the charge 
measurement is very high, the phase becomes completely 
undetermined. The information that can be gathered 
from the SET contains no information of the phase and 
thus the precision of measuring the charge within the 
qubit (|a| 2 and |6| 2 ) can be very high. Furthermore, 
since charge is not disturbed by the measurement for 
imix tmeas; w e are essentially doing the classic quan- 
tum non-demolition measurement. 

The time scales are calculated for a specific set of pa- 
rameters. With aluminium structures in mind, the su- 
perconducting gap of the SET is chosen to be A = 2.3 K, 
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Figure 3: The probability P(m,t) that m electrons have tun- 
nelled in a superconducting SET plotted as a function of time 
t for initial amplitudes V0.75 (n — 0) and \/0.25 (n = 1). The 
parameters are for aluminium structures and they are given 
in the text. The graph shows the separation of the two peaks, 
the faster corresponding to the qubit's state and the slower 
to the state 1. The two peaks are clearly distinguishable from 
each other, and thus the measurement time is very small. 
The curve in the box shows the corresponding evolution for a 
normal-metal SET. 



time and no separation of peaks is visible. Thus, for alu- 
minium structures, the normal-metal SET is unlikely to 
be sensitive enough quantum measurement device in the 
high-voltage regime. Both for the normal-metal and su- 
perconducting SET the qubit state |1) or |0) determines 
the current through the SET by shifting the energy lev- 
els in the middle island. However, due to the divergence 
in the superconducting density of states at the gap this 
effect is strongly enhanced and leads to almost complete 
blocking of the current for n = 1 in the superconduct- 
ing case. In general, this indicates that properties of 
the superconducting state can be very useful in perform- 
ing sensitive quantum measurements. In particular, the 
superconducting SET can be considered to constitute a 
very good quantum nondemolition measuring device as 
the mixing time is much larger than the measurement 
time imixAms 3> 1< for parameter ranges accessible by 
aluminium structures. 
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the capacitances Cj = 8.0 • 10- 16 F, C T = 8.0 • 10- 17 F 
and all the remaining capacitances are set to 1.0 • 10~ 17 F. 
These fix the charging energies E q ^ w 1.0 K, _E sct w 5.2 K 
and E[ nt w 0.13 K. The gate charges and tunnelling 
coefficients are chosen as Q q b = 0.35, Q sot = 0.15, 

\T L \ 2 = \T R \ 2 = 0.0025. To justify the charge-qubit 
approximation E q ^ ^> Ej the strength of the Josephson 
coupling is chosen as Ej — 0.05 K. Finally, the bias- 
ing voltage is chosen as eV = 16.5 K. This set of val- 
ues gives the measurement time t ms « 3.6 • 10~ 10 s, the 
dephasing time t$ ~ 7.9 • 10~ 10 s, and the mixing time 
imix ~ 2.3 • 10 _7 s. The ratio between the mixing and 
measurement times is very high, i.e. i m ix/ims rj 625. 
The same set of parameters for a normal-metal SET 
(A = 0.0 K) gives t ms « 2.6 • 10" 6 s, « 9.32 • 10- 10 s 
and i mix « 1.0 • 10~ 7 s. 

The system has been simulated numerically for the pa- 
rameter values above, and the probabilities P{m, t) for m 
electrons having been tunnelled during time t are plot- 
ted for the superconducting and normal-metal cases in 
Fig. | In the case of the superconducting SET, the mea- 
surement time is very small, as only one peak (corre- 
sponding to the qubit's state n = 0) propagates in time 
and the other peak (n = 1) is blocked by the presence of 
the gap. The mixing effects are small, as the peaks re- 
main separated and the small visible spreading is caused 
mainly by the shot noise. For the normal-metal SET 
the measurement time is much longer than the mixing 
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The Basset function of integer order K n is gener- 
ally defined as an elliptic integral [[l5| K v (z) = 
Ji°° i 2 tu/2-u ■ The numerical analysis of 



r(i+„) \2) Ji (3,2-1)1/2- 
the Basset function for small values of the argument 
z is based on presenting it as a Neumann series of 
Bessel functions. For a large argument, a good al- 
gorithm can be written by applying the Pade tech- 
nique on the asymptotic approximation for large z 

K v (z) ~ V^-[l + Sr=i ( - 1/2 ^ +1/2) - ], where 



x(x + l)(a; + 2)- 



(x + j — 1) is the Pochhammer 



polynomial. 



